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Abstract 

We study the statistical mechanics and dynamics of crystalline films with 

a fixed internal connectivity on a random substrate. Defect free triangular 

lattices exhibit a sharp transition to a low temperature glassy phase with 

anomalous phonon fluctuations and a nonlinear force-displacement law with 

a continuously variable exponent, similar to the vortex glass phase of directed 

lines in 1+1 dimensions. The periodicity of the tethered monolayer acts like a 

filter which amplifies particular Fourier components of the disorder. However, 

the absence of annealed topological defects like dislocations is crucial: the 

transition is destroyed when the constraint of fixed connectivity is relaxed 

and dislocations are allowed to proliferate. 
46.30Pa 68.35Rh 76.60Ge 
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I. INTRODUCTION 



The pinning of elastic media subjected to external forces is a subject of considerable 
interest in connection with a variety of phenomena, ranging from epitaxial surface growth to 
transport in planar arrays of flux lines in type II superconductors |l|]. High temperature su- 
perconductors are especially interesting in this respect because of the simultaneous presence 
of large thermal fluctuations and of quenched disorder. Much effort has been devoted to the 
study of 1+1 dimensional models, which are models of vortex lines confined to a plane 
(see Fig. la). Although some quantitative questions have yet to be answered satisfactorily, 
a clear qualitative picture of the physics involved has emerged, which can be summarized 
as follows: the 1+1 dimensional flux array, subjected to external point disorder, displays 
a transition, at some temperature Tg, between a high temperature regime, dominated by 
thermal fluctuations, and a low temperature regime, where the behavior of the system is 
controlled by a line of fixed points. This is a transition to a disorder-dominated phase, where 
the elastic system is pinned. Pinning affects both static correlations and dynamic responses 
in a nontrivial way, giving rise to nonlinear current-voltage characteristics. Crucial to pin- 
ning is the discrete nature of the elastic system, which, roughly speaking, acts as a Fourier 
filter for components of disorder on length scales corresponding to the lattice spacing (i.e., 
the distance between flux lines). 

The question naturally arises, whether coherent amplification of the disorder (assumed 
to exist at all physically relevant wave lengths) is instrumental to the pinning of a vortex 
array; that is, it is important to question the role of long range crystalline order of the array 
in selecting out particular Fourier components of the pinning potential. This feature is built 
into the 1+1 dimensional vortex line model, which possesses algebraic crystalline order at all 
nonzero temperatures and is also topologically perfect. The topological perfection arises for 
vortex lines because the average magnetic field is parallel to the plane, and thus the lines 
cannot terminate in the plane, and their labelling in a perfect crystalline sequence is always 
unambiguous. 
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Alternative two dimensional models of elastic solids can be considered, which allow for 
topological defects, such as dislocations. These are 2+0 dimensional models (see Fig. lb). 
For vortices in superconductors, the average magnetic field must then be perpendicular 
to the plane of the film, and defects leading to multivalued displacements of the vortices 
are allowed, as they are in many other experimental situations. Additional experimental 
realizations include colloidal crystals |^ , amphiphillic mono- or bi-layers composed of lipid 
molecules (possibly polymerized) p, electrons in semiconductor heterostructures P], and 
magnetic bubble arrays [0]. The constraint of fixed nearest neighbor connectivity could be 
enforced in some cases by polymerization or, more generally, simply by large kinetic barriers 
to particle exchanges at low temperature. These systems can be topologically perfect, or else 
can be subject to either quenched or annealed internal defects, in addition to the external 
pinning potential. 

The goal of the present work is to develop the theory of pinned two dimensional crystals 
within the framework of two dimensional elasticity theory. For the case of a topologically 
perfect crystal of identical tethered particles, subjected to thermal fluctuations and quenched 
pinning, the behavior is qualitatively similar to the 1+1 dimensional model (displaying, e.g., 
a sharp phase transition to a low temperature pinned phase). The possibility of topological 
defects, however, introduces significant modifications. 

Annealed dislocations destroy positional quasi long range order in two dimensional crys- 
tals above some finite melting temperature Tm 0, where a liquid crystalline "hexatic" phase 
exists, with algebraic long range order in the bond angle. As will be shown below, the melt- 
ing temperature is always smaller than Tg, so that a two dimensional crystal is always melted 
at the temperature below which pinning disorder would become relevant for a topologically 
perfect solid (i.e., dislocations are a relevant perturbation at Tg). Thus, the phase transition 
discussed above is washed out in the presence of thermally excited dislocation pairs. Anal- 
ogous behavior is encountered in the random field XY model 0]. An extension to vector 
displacement fields has been studied by Giamarchi and LeDoussal [Q. 

The hexatic liquid crystalline phase of the untethered membrane above Tm (see Fig. 2a) 

3 



displays similar behavior when subjected to a component of the random substrate disorder 
which couples directly to the bond angle field. The analogy with the random field XY 
model becomes a rigorous mapping for annealed hexatic membranes: either disclination 
unbinding or substrate disorder always distabilize the hexatic line of fixed points, and it 
is unclear if there is a sharp finite temperature phase transition. Polymerized tethered 
membranes behave differently, however. Although quenched-in unpaired dislocations destroy 
translational long range order, they cannot drive the shear modulus to zero. The finite shear 
modulus makes the bond angle fluctuations "massive" ^j. These fluctuations are now stable 
to weak external disorder. 

Quenched topological disorder has recently been studied for "tethered surfaces" |T0| , 



in which defects are frozen into a two-dimensional network of covalently bonded particles 
fluctuating in three dimensions. It is of considerable interest to determine what happens 
when such disordered tethered surfaces are forced to lie flat and brought into contact with 
a disordered polycrystalline or amorphous substrate. A particularly simple example of such 
tethered disorder is shown in Fig. 2b, where a topologically perfect triangular lattice is 
disrupted by random substitutional disorder. Unfortunately, the method used in this paper 
cannot be directly applied to such systems ||TT| . Tethered substitutional disorder invalidates 



a straightforward analogy with random fleld models (see Sec. 2). Cule and Hwa have studied 
this problem in one dimension, and concluded that a new, stronly pinned glassy phase arises, 
characterized by exponents in the "random manifold" universality class [0 . Similar results 
may apply to two dimensional tethered networks with quenched-in vacancies, interstitials, 
dislocations or disclinations as well as random substitutional disorder |12| . 



A theory of three dimensional tethered networks with both quenched random internal 
defects and a quenched random external potential would have interesting implications for 
the tangled arrays of vortex lines which may arise when bulk Type II superconductors are 
subjected to strong external magnetic flelds. If melted flux liquids are cooled rapidly, barriers 
to flux cutting |TB[ may become sufficiently large that the vortex lines freeze at low temper- 



atures into a nonequilibrium directed "polymer glass" [|lj] . The usual triangular Abrikosov 



flux lattice would then be disrupted by a quenched array of dislocation and disclination 
lines of arbitrarily large size. The resulting vortex array could have a shear modulus over a 
wide range of experimental time scales because of entanglement constraints, and would be 
subject to point pinning by imperfections in the underlying host superconducting material. 

This paper is organized as follows. In Sec. II, we present the theory of ideal (topologically 
perfect) two dimensional crystals subject to external point disorder. Effects of annealed dis- 
locations are discussed in Sec. Ill, and some concluding remarks are presented in Sec. IV. 
Some of the effects of quenched substitutional disorder are discussed in Appendix A. Tech- 
nical details of the derivation of the renormalization group recursion relations used in Sec. II 
are contained in Appendix B. Appendix C derives the exponent z for a relaxational model of 
tethered crystalline membrane dynamics. An Ornstein-Zernicke description of the hexatic 
phase is derived in Appendix D. 



II. PINNING OF IDEAL CRYSTALS 

Consider a two dimensional Bravais lattice Rmn = + nb and denote by u{m, n) the 
displacement from the equilibrium position Rmn- We begin by restricting our analysis to 
an ideal lattice of identical particles, where the displacements are single- valued functions of 
position. Thus, we exclude for the moment the possibility of topological defects. The strains 
associated with random substitutional disorder will be discussed later in this section. The 
energy of an ideal crystal undergoing a small deformation can be expressed, using continuum 
elasticity theory, by the harmonic hamiltonian 



H, = ]^j (fr{\ul + 2iiul), (2.1) 



where A and /i are Lame coefficients and the strain tensor is defined as 



"ii(^) = nid,iUj + djUi). (2.2) 



The only regular lattice in two dimensions with sufficient symmetry to be described by the 
isotropic elastic theory ( p.l|) is a triangular array p|. In this case, it is convenient to take 



a = ao(l,0), ^=^oQ,^), (2.3) 

where is the lattice constant. 

If the lattice is subjected to an external pinning potential, the total pinning energy is 
given by Upin = J2mn V[Rmn + u{m, n)]. With the aid of Poisson's summation formula, we 
write 

Upin= 1 d^rV{r) ^5'^[f- Rmn - u{m, n)] 

mn 

= I (frV{r) I dai f da^S^ir - a^a - - u{ai,a2)]Y.e^'''^P'''+''''^\ (2.4) 

where {m,n) and {p,q) are pairs of integers. Next, we change integration variables from 
cTi, (72 to r = (x, y) = aid + + u{ai, a2). The Jacobian of this tranformation is 

daida2 daida2 1 , - , . 

dx dy dy dx \a x h\ 

The 5 function in Eq. ( P^ ) fixes the value of a to be the root of 

aid + 02h + u{ai,a2) = r, (2.6) 

which can be expanded in series of the small displacement u. Upon noting that the zeroth 
order term is 

^2 = — ^ > (2.7) 



\d X 6| |a X h\ 

and using the 5 function to eliminate the integrals over ai and (72, we can write the total 
pinning energy as 

Upin ~ / d^rV{r)^-^{l -V-U)Y^ ^iG^^-[r~u{r)] ^ (2.8) 

J Ia X /)l ^„ 



mn 



\d X b\ 

where the {Gmn} are reciprocal lattice vectors. The total hamiltonian of the system is 
obtained by adding the hamiltonian of an ideal two dimensional crystal to the pinning 
energy. We study the effect of a random distribution of weak pinning potentials, with mean 
and variance defined by 



V{f) = 0, V{r)V{P) = A5^{f- r'). (2.9) 

If we restrict our attention to Gmn = and to the six smallest reciprocal lattice vectors 
in the summation of Eq. (5), then the total Hamiltonian which comprises our model may 
be approximated (up to an additive constant) by 

n=yd\^Xul + 2fiul-w{r)uu + j2VG,{f^e-''^-''^'^y (2.10) 

where the Gi, Z = 1, 2, 3, are three reciprocal lattice vectors inclined at 120° angles to each 
other in the innermost ring, 

w{r) = XllX, (2.11) 
\a X b\ 

and V^{r) is a local Fourier component of the random potential 

The integration is carried out over an area AQ, centered on r, large compared to the lattice 
spacing but small compared to the sample dimensions. We have neglected terms of the 
form V ■ uexp{iGi ■ [r — u{r)]), which are less important than the ones we have kept. In 
Appendix A, we show that random substitutional disorder contributes to w{f) and leads as 
well to a random term of the form — ^ / (Prwij{r)uij{r) . A two dimensional bead and spring 
model with random spring lengths contains similar contributions. Quenched substitutional 
disorder has important additional effects, however. Internal disorder in the particle sizes or 
bond lengths violates the discrete translational invariance of the Hamiltonian (2.10) under 
u{r) — s> u{r) + Rmn- Cule and Hwa have argued that this symmetry-breaking leads to an 
effective "random manifold" potential which depends in a complicated way on both f and 
u{r) , and causes an instability to a more strongly pinned glassy phase than the one studied 



here \12 



For tethered networks of identical particles, we expect that pinning effects due to the 
random phases and amplitudes embodied in the {Vg,} will be important below some critical 
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temperature, and that the properties of the system in the vicinity of this temperature will 
be perturbatively accessible by renormalization group methods. 

The starting point of the renormalization program is perturbation theory. Computation 
of the disorder- averaged observables of interest, such as the free energy or the two-point 
functions, requires expanding the logarithm of the partition function Z in powers of the 
weak pinning potential and averaging term by term. Such averages are conveniently handled 
by the rephca trick, which involves first calculating = l-l-nlogZ-l-O(n^) and eventually 
taking the limit n — > 0. The disorder average of leads to 

Jvui...Vune-^-^, (2.13) 

where the harmonic term is 

-Y = ^EJ j^,u,^i-k)k\[^iP^,+{2^i + X)P^^^^^^^ (2.14) 
and the interaction term is 

^ = -^E / d'r±cos[G, ■ {u^{r) - (2-15) 

1=1 •' a/3 

Greek indices label different rephcas, while latin subscripts = 1,2) are used for the com- 
ponents of the displacement vector u. Summation over repeated indices i,j is understood. 
The vectors Gi are the smallest nonzero vectors in the reciprocal lattice. For a triangular 
lattice of spacing ao, / = 1, 2, 3 and \Gi\'^ — 167r^/3aQ. Effects due to reciprocal lattice vectors 
of larger norm are irrelevant. 

The structure of the replicated Hamiltonian is quite simple. The harmonic part contains 
a term diagonal in the replica indices. This is simply the replicated hamiltonian of an ideal 
two dimensional crystal in Fourier space, where and P^ are longitudinal and transverse 
projectors, respectively. In addition, transverse and longitudinal terms are considered, which 
are constant in replica space. The tranverse term, while not present initially, is generated 
by renormalization. Its coefficient B can be set to zero initially. The cosine term arises as a 
consequence of the discrete nature of the lattice. Its amplitude is related to the correlation 
function of the randomness by 
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where f2 = |a x 6| is the area of the unit celL 

The use of rephcas brings out an important property of the model. Note that the 
interaction term involves only differences of fields with different replica indices. Hence, the 
"center of mass," in replica space, of the fields m^, 

^ = ^-1/2^^^, (2.17) 

a 

is a free field and does not suffer renormalization. This symmetry implies that 2/i + A — 
nA, fi — nB do not renormalize, a result valid to all orders in perturbation theory, similar 
to the invariance under renormalization of the spin wave stiffness in the random field XY 



model llg. As a consequence, the renormalization group flow of the disorder coupling, is 



one dimensional. 



Next, consider the connected Green's function 



< Ui{r)uj{f') >, 



dJi (r ) dJj (P) n dJi (r ) dJj (F 
obtained from the generating functional 



n=0 



(2.18) 



= < exp(^i j d?rJ{r) ■Y.^c.ir)) >R ■ (2-19) 

Hereafter, the notation < >r will stand for average with respect to the integrand in 
Eq. ( p.l3| ). Since the source J couples only to ^, a free field, this Green's function is 
the free correlation function, independent of g, 

< Ui{r)Uj{r') >c = -Sij- — - — — - log h const, |f - fo| oo, (2.20) 

47r/i 2/i + A ao 

where the connected part is < A{x)B{x') >c=< A{x)B{x') > — < A{x) >< B{x') >. 

The peculiar properties of the glassy phase do appear, however, in the nontrivial behavior 



of some response functions as well as in the full correlation function < Ui{r)uj{f') >, which, 
according to the analysis above, is a probe of sample-to-sample fluctuations. In the language 



of replicas, these fluctuations are captured by introducing a replica dependent source field, 

— * 

Ja{^i and a new generating functional 

Z-[{JJ] =Z^<exp(z j (fr^J^{r)-u^{r)^>R. (2.21) 

— * 

Differentiation with respect to Ja yields the Green's functions 

Gijap{r- f^) =< Uio,{r)Uj^{7^) >r . (2.22) 
Provided symmetry under permutation of the replica indices holds, we can write 



< Ui{r)Uj{r') >c= lim[G'ijii(f - r') - Giji2{f- f)] 



< Uiirjujif) >= lim G',,n(r - f), (2.23) 

where the limit n ^ is here simply a convenient bookkeeping device for doing perturbation 
theory. 

The perturbation series for the disorder-averaged Green's functions diverges in the ther- 
modynamic limit at low temperature. The divergent diagrams are most easily recognized 
by considering the expansion of the free energy, 

= - lim i (- < Hi >K +^(< {nir >n-< (Hi) >%) + ..). (2.24) 

Upon defining a reduced temperature 

= (2.25) 
one finds, up to regular terms, in order g: 

_ < Hj >R 

T 

and in order (g^): 



gca^n{n - 1)?>t: {L / ay/cf^ , (2.26) 



2^ < {Hif >u- -g\^a'n{n - l)[IoW + {n- 2)/o(V'/2)]37r^(^)'^^((L/aV^)^- - 1 

(2.27) 
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Infrared and ultraviolet cutoffs L and a, respectively, have been introduced; c = \e^^ ~ 0.79, 

T\G\\^i+\) 

47r/i(2/i+A) 



where E is Euler's constant; Jq is a modified Bessel function, and ip = V'^iy'l^) ■ '^^^ details 



of the calculation can be found in Appendix B. 

The divergences can be removed order by order in a double expansion in powers of g 
and T. The parameters of the renormalized theory transform, under rescaling of length by 
e\ according to the following equations: 

dX 

I- 

^=2r~g-^f[2Ioi^/2)-Iom 

where g = ?>Txgcc?{ljl a^fc)^'^ . The flow of the disorder coupling to zero for r < 0, i.e., for 
T > Tg, means that the discreteness of the lattice is irrelevant in the high temperature 
phase. Thus, the correlation functions in this phase are similar to the gaussian model of 



Ref. [1171. particular. 



< Ui{r)Ui{0) > ~ -r^logr (2.29) 

where 

_ T 3fx + X A A B \ 

Below Tg, on the other hand, the disorder coupling flows toward a finite fixed point value 
of g* = 3t/[2Iq{iIj/2) — Iq{iIj)]. The runaway flows of A and B cause the correlation function 
to grow as log^ r, a behavior which was termed "superroughening" in the context of scalar 
models of surface growth |]I8[ . 

As we anticipated, another distinctive property of the glassy phase is the nontrivial near- 
equilibrium dynamics. The dissipative dynamics of the system embodied in the Langevin 
equation 
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with ( a thermal noise, can also be studied by dynamical renormalization group methods 
||T9| . The detailed treatment of this model is described in Appendix C. Regularization of the 
perturbative expansion of the dynamic response leads to a renormalized friction coefficient 
7, from which the dynamic exponent can be extracted 

z=2, T> Tg 

24 2/i + A / /i \3jiT\ T 
'^'^^V^c3j^[2j^) [2/o(^/2)-Jo(^)]' ^ 

where c = |e^^ ~ 0.79 is the same function of Euler's constant as appears in the static 
calculation. 



III. EXTERNAL DISORDER AND TOPOLOGICAL DEFECTS 

The foregoing discussion assumed the fixed topology of an ideal lattice. Spontaneous 
nucleation of topological defects, which occurs above the melting temperature 

destroys translational long range order 0]. The ratio of the glass temperature to the melting 
temperature is always greater than one (in fact, T^/Tm > 6 for all elastic constant values 
in the physically relevant range /i > 0, /i + A > 0), so that dislocations are expected to 
be a strongly relevant perturbation at Tg. Thus, either the random substrate potential or 
thermally excited dislocation pairs always destabilize the harmonic hamiltonian, and the 
transition of the previous section does not occur in the presence of annealed topological 
defects. 

Although annealed dislocations destroy translational long range order, the resulting hex- 
atic phase does possess (algebraic) long range order in the bond angle. A harmonic hamil- 
tonian for the hexatic phase can be obtained in Ornstein-Zernicke approximation, valid at 
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long wavelength. The details of the derivation are presented in Appendix D. We can con- 
sider the stability of the long wave length hamiltonian, Eq. (D14), to an external random 
potential coupled to the bond angle. An experimental realization of this system is pro- 
vided by a hexatic liquid crystalline film adsorbed onto a polycrystalline substrate, that is, 
a substrate whose randomly varying crystallographic axes locally bias the orientation of the 
bonds in the film. The total hamiltonian becomes precisely that of a random field XY model, 
which was studied by Cardy and Ostlund P]. Similar to the case of untethered crystalline 
films discussed above, it follows from ref. that the harmonic hamiltonian (D14) is always 
destabilized either by the external disorder or by thermally excited disclinations. There is 
an important difference between crystalline and hexatic membranes, however. In the case 
of a crystalline membrane, the ideal topology can be fixed by polymerization, and the fixed 
line discovered in Sec. II should be experimentally accessible. In contrast, it is impossible to 
prevent disclinations from destabilizing the vortex glass fixed line in a hexatic film, except 
by quenching the topology of the film. But in this case, Eq. (D14) reveals, upon treating the 
singular density and bond angle fluctuations, 5ps and 56 as quenched variables, that this 
quenched hexatic phase has a finite shear modulus, rendering the bond fluctuations massive. 
Because of the finite shear modulus, the bond angle excitations are stable to weak disorder. 
Whether a tethered hexatic (or liquid) film is unstable to a "random manifold" glassy state 



|12| is an interesting subject for future investigation. 



IV. DISCUSSION AND CONCLUSIONS 

In this paper, we have studied the physics of 2+0 dimensional arrays of identical particles 
in an external random potential, using a vector extension of the 1 + 1 dimensional random 
phase model. The predictions for vector models without topological defects are qualitatively 
similar to the scalar model: the discreteness of the vortex array coherently enhances the 
Fourier components of the external disorder which are commensurate with the lattice, leading 
to a glassy phase at low temperature. This phase is characterized by a nonlinear response to 
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an external driving force, as well as by static correlations which diverge more strongly than 
a simple logarithm. Quantitative expressions for the static and dynamic exponents near the 
transition were computed within a perturbative renormalization group scheme. 

The physics changes completely, if annealed topological defects are allowed in the two 
dimensional lattice. The possibility of important topological defects constitutes the prin- 
cipal difference between the 2+0 and the 1+1 dimensional models of vortex arrays. Both 
annealed and quenched dislocations have been considered, the most interesting case being 
provided by quenched dislocations. The nonvanishing shear modulus of a membrane with 
quenched-in dislocations prevent the bond angle from following the random bias of an exter- 
nal polycrystalline substrate, so that the bond angle order parameter in quenched hexatics 
is stable to weak disorder of this type. 

Our study is relevant to several other situations besides the pinning of vortex arrays in 
type-II superconductors. Systems of current experimental interest were mentioned in the 
Introduction. Here, we would like to comment on possible applications to tribology, the 
study of friction and lubrication. We are interested in the behavior of two surfaces brought 
in contact and rubbed against each other in the presence of an intermediate thin layer of 
lubricant. This boundary layer is often modelled as a two dimensional, incommensurate 
crystalline overlayer ||2l]|. Our work may be useful in generating more realistic descriptions 



which allow for a) surface imperfections, acting as pinning centers on the lubricant overlayer; 
b) changes in topology of the overlayer, especially excitation of dislocations, which must 
surely be important at finite temperature and/or under finite stresses. We leave the pursuit 
of this interesting topic to future work. 

Noted added: After this paper was submitted, we received an interesting preprint by D. 
Carpentier and P. Le Doussal ( |cond-mat / 96 1 1 1"68| ) which reaches similar conclusions using 



a different renormalization group method. Comparison with their results enabled us to 
uncover an error in the first version of our paper which, although it did not affect our basic 
conclusions, changes the coefficients in our recursion relations. Once the error is corrected, 
results obtained by the different methods agree. We are grateful to P. Le Doussal for bringing 
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this discrepancy to our attention. 
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APPENDIX: A 

Effects due to a disordered substrate were incorporated into isotropic two dimensional 
elasticity theory in Section II. In this Appendix, we discuss effects on the elastic properties 
due to random substitutional disorder in the membrane, as exemplified by the large impurity 
atom displayed in Fig. 2b. We do not discuss the important interplay between random 
substitutional disorder and the disorder substrate potential |T^. A distribution of impurity 
atoms with sizes different from the average leads to random strains. Annealed defects of 
this kind can be integrated out and simply alter the elastic constants fi and A. As we shall 
see, the strains in the quenched case contribute to the coefficients A and B displayed in the 
replicated hamiltonian, Eq. ( p.l4|) . Quenched random vacancy or interstitial defects as well 
as tightly bound dislocations pairs or triplets would affect A and B similarly. We work with 
a continuum model studied already in the context of random tethered surfaces fluctuating 



in three dimensions |jTO[ . The only change required is neglect of displacements normal to the 
average plane of the membrane. These phonon modes become massive due to the interaction 
with the substrate and can be integrated out without affecting our basic results. 

We assume a topolgically perfect lattice and replace Rmn by a coarse-grained function 
R{f) which gives the lattice displacement i? as a function of the reference position r. We 
use a generalization of Eq. (|2.1|) , 

H=\j d'r{Xul + 2^^u%), (Al) 
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where the strain matrix is now given by [10 

= ^ {diR -djR- diR° ■ djR^). (A2) 

Here, R^{r) is a preferred lattice distortion which minimizes the energy in the absence 
of thermal fluctuations. In the absence of defects, diR^ ■ djR^ = 6ij. Localized defects like 
substitutional disorder, vacancies, interstitials, etc., lead to deviations which we parameterize 
by 

diR"" ■djR'' = 6ij + Ci,{f). (A3) 
If we assume uncorrelated gaussian disorder, the probability distribution of the tensor Cij{r) 



takes the form 10 



Vr[ci,ir)] oc exp(-^ / d^cl - (frcl^. (A4) 

We now set R{r) = r + u(r), so that 

diR ■ djR = 6ij + ^{diUj + djUi) + ^diU ■ djU ^ 6ij + Uij. (A5) 



The hamiltonian ( |Al| ) then takes the form 



H = const. + 2 y '^^^(•^'"ii + '^f^'^fj ~ •^Cji(r)'Uj.,(r) — 2fiCij{f)uij). (A6) 

The term proportional to Uu represents random dilations or contractions due to isolated 
impurities in positions of high symmetry, while the more complicated tensorial coupling 
describes more anisotropic defect configurations. Upon replicating this hamiltonian and 
tracing out the Gaussian disorder, we obtain contributions to the coefficients A and B in 
Eq. (ra . 



APPENDIX: B 



This Appendix details the calculations leading to the flow equations for g, A, and B. We 
begin by evaluating the right hand side of Eq. (|2.24| ) term by term. To order g, 

16 



exp [—11- 

Before proceeding with the calculation, we must introduce cutoffs to deal with infrared and 
ultraviolet singularities arising from integrals of the type 

I J2^^''^^"^' ^^^^ 

A long wave length cutoff L is introduced to eliminate infrared divergences. Its effect 
amounts to shifting l/fc^ — > + L^^). The limit L ^ oo can be taken safely at the 

end of the calculations. The ultraviolet divergences are removed in coordinate space by the 
simple shift r — > \^r^ + a^, where a is a short wave length cutoff of order the lattice constant. 
Note the short distance limit: 



2u „ik-s 



(27r)2 P + L- 



1 ^ .Vf^T^ 



1 (r^ + a^) r , 



Here, Kq is a modified Bessel function, and c = (l/4)e^''', where 7 is Euler's constant. Thus, 
the asymptotic behavior of the integrals in Eq. (|B2D is readily evaluated 



(27r)2 * A;2 ^ Sir \ Anr- stt 

Upon substituting into Eq. (pi]), and recalling Eq. (|2.25|) , which defines the reduced tem- 



perature r = 1 - J = 1 - ^fgx' ^6 obtain Eq. ( g:26D 



Next, in order g"^, one has 

\ 2\ 3 3 n n 



-'■ ^ ' 1=1 l'=l afi a' 13' J •' \ ^ 

xexpfe-«(0-«;3(0) 



3 3 n »i / \ (4-4r) 

^^eeee(-: 
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Svr^ 2/i + A L2 



where = 5q,q,/ + 6ppi — dap' — ^a'p, and use has been made of Eq. ( p3D in anticipation 

of the thermodynamic hmit. Hereafter, a = y/ca and ip = 

— * 

Logarithmic divergencies appear, at the critical temperature (r = 0), when Gi ■ 
^I'Qap /I^P = 1- The logarithmically divergent terms contribute to the multiplicative 
renormalization of the coupling g, while all other terms can be absorbed into an additive 
constant. (Had we chosen to renormalize the self energy, the need of an additive renor- 
malization would not have arisen, but the calculations would have been considerably more 
involved.) To proceed, we first fix /, a, and (3, which can be done in 3n{n — 1) indepen- 
dent ways. Then, logarithmic contributions arise either for /' = /, Q'^p = 1, or for I' ^ /, 
= —2. The respective combinatorial weights are 2(n — 2) and 2. Carrying out the 
remaining integrals leads to Eq. ( |2.27| ). 

We now compute the flow of the coupling constants A, B under renormalization group 
transormations. Because the self energy is momentum independent in the first order of 
perturbation theory, the lowest order contribution to the perturbative renormalization of A, 
B arises in 0{g^). The free propagator in momentum space, obtained from Eq. ( |2.14| ), can 
be written as the sum of two terms: 

i\T ^ ^ \ ^ - rRfil 

' ^^{^^-nB) '^^ k^{2^l + \){2^l + \-nA) ''^ ^ ' 

The perturbative expansion of the generating functional Eq. (|2.21|) through second order 

reads as 

Z-[{Jj]=exp(-| j -fLY^Y.j^^^k)M,,ik)Jj^i-k) + Y.J,^^^^ 



a/3 

2 



x[l + gJ^,{J) + ^J^,{J) + ...] (B7) 



where 



X exp 



flA(4-4r) 



"' Q/3 a'13' '' ^-^^ ^ J 



Thus, the renormahzed propagator N^{k) can be calculated perturbatively as 

(27r)V <^'^2(J) 



TNm)S{k + k') = -TNij{k)6{k + k') + 



2 5Ji^{k)SJj^{k') 



+ 



with 



5Ja(A;)5J,>(A;') 



><7fuIl(G!'M,i(k)(5^^-5f,^)e-''^U^^ 



x( G'^,Mj,j(k')(5^, - 5p,)e~"''-''+ G''^,My^{k'){8^,, - 5p,,)e 



Contributions behaving as l/Zc^ for small k arise from 

EEE p^r^v(f )^""exp(-Tg^|'EG^^.(r - 0^5' 



j'2 



5\(4-4r) 



M-TQfpY.G\M,,{s)G^; 



^3 



_rp2 



^ mm' 



(B8) 



(B9) 



(BIO) 



(Bll) 



Power counting reveals that logarithmic divergencies arise when I — I' and a — (5' and 
(5 — a' . Thus, with the help of the following relations. 



E^m^m' - -^\G\ ^mvnJ 
I ^ 

i '~ " 2 3 2 T L 



(B12) 
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we obtain 



3n 



2 V/x2"^^^[^°^^^ + 2^^^^^] 



+ 



1 



(2/i + A)2fc2 



4r 



4r 



1 



(B13) 



Upon projecting out of Eq. ( [B13| ) the longitudinal and transverse components, and recalling 
that IGpT = 87r/i|^ji| + 0(r), we finally arrive at the fiow equations 



4r 



rf5 



4t 



(B14) 



APPENDIX: C 

This Appendix details the calculation of the dynamical exponent z. Consider the 
Langevin equation for the free hamiltonian Hq of Eq. (|2.14| ) 

where C is a white noise defined by the correlation 

< >= 2D5,,5{r- 7^)5(1 - t'). (C2) 

The fiuctuation-dissipation theorem insures that thermal equilibrium at temperature T is 
estabhshed for D = T7. 

Equation (|C1D is easily solved. The momentum and frequency dependent transverse and 



longitudinal correlation functions read 



The presence of an external pinning potential turns Eq. ( |CI| ) into Eq. ( |2.31| ) by intro- 
ducing on the right hand side the additional term 
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which can be regarded as a perturbation. First order perturbation theory can thus be 
employed to compute the correlation functions. The perturbation series is found again to 
diverge logarithmically below T^. This divergence causes the dynamical exponent to deviate 



from its mean field value, z = 2 \22\. To quantify this deviation, consider the transverse 
correlation function 

1 



< Ui{k,uj)uj{—k, —uj) >^ 



+ Pl,Pl< Fmik,Uj)Fni~k, -iu)>] + 



(C5) 



where the overline indicates averaging over the external pinning disorder and we have defined 



Jujt-ik-f ^Upin 



6Ui{f,t)' 



Upon carrying out the disorder average, we obtain 

3 



<Fi{k,uj)Fj{-k,-uj)> = fdte'''%Y.G\G^j <cosG'- (M(f,t) -M(f,0))> . 

■' 1=1 

Noting that 



(C6) 



(C7) 



a^2-2r 



<cos&-{u{0,t)-u{0,0))>={j) exp(x:/|;iG^(5°^(^,t) + ^S'(^'^))^5 



a^2-2r 



exp 



\G\^T roc 



_2£i±Afc2|j| 

e 'I 



/e T 

\ ^ik^ (2/i + A)A;2 



d ^ 2-2r 



L2 



and setting = 0, we arrive at 



<Fi{k,uj)Fj{-k,-uo)> = 2gd'—\GY6., 



2/x' 



2/i + Ay \2^/c'^-^ ^i\t 
/i 27 



l-T 



(a 



2/x + Ay 3/i + A a/c 



L 



log - +0((7r). (C9) 



Thus, we can write 



L-— = 2A'Kga _ 
ctL Va/ 3/i + A\2/i + A. 

The dynamical exponent is determined by the value of C,{g) at the fixed point 



(CIO) 



(Cll) 



leading to Eq. ( ^ ). 
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APPENDIX: D 



The Ornstein-Zernicke theory of hexatics is derived in this Appendix. We start from the 
partition function of the crystal 

Z = Jv[u,]V[u2] exp(-F/T), (Dl) 

where 

F=yd'r{Xul + 2^^ul). (D2) 

In the presence of dislocations, the displacement field u is multivalued; if an arbitrary loop 
C encloses dislocations with total Burger's vector b, then 

diUjdxi = bj. (D3) 

It is convenient to express Eq. ( |D3| ) in local form, 

e^kejl^k^lu,,{f) = ^ x V5(r - f„) = S{r}. (D4) 

Here, we have introduced a source term, S{r), related to the local density of Burger's vectors; 
it can be easily generalized to include other types of defects, such as disclinations, vacancies. 
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The partition function can now be evaluated as an unresticted integral over the three 
independent components of the symmetric strain tensor Uij, provided that the constraint 
is enforced. This is easily accomplished with the aid of an auxiliary field ip: 



Z = Jv[un]V[uu]V[u22]V[tlj]exp{-F/T)exp(i J d^r^{r){eikejidkdiu,,{r} - S{f))y (D5) 

Next we write, as can be most generally done for any symmetric tensor, 

u.,{r-) = ^id.v,+d,v,) + P;'^h, (D6) 

where v{f) is a single-valued vector (describing, e.g., the small thermal displacements of 
the atoms from the sites of a randomly distorted lattice), and h{r) a scalar field, related to 
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the defect density S{r), in terms of which the constraint enforced by the auxihary field i/) 
becomes 

V^h{r)=S{r). (D7) 

The partition function becomes a functional integral over v, h, and ip, and the latter 
fields can be integrated out easily to obtain the free energy functional 

+A J d\{V ■ v)V~^S{r) + + A) 1 dh{V-^S{r))\ (D8) 

This is a rather general expression of the two dimensional elastic free energy. It allows us 
to study the effect of both quenched and annealed defects. 

First, we consider the effect of annealed dislocations. It is instructive to recast Eq. (P8|) 
in terms of density fluctuation and bond angle variables, 5p and 9 respectively, which are 
related to the divergence and curl of the displacement vector: 6p = V ■ u, 6 = V x u/2. 
Because the displacement vector decomposes naturally into single- valued and singular parts, 

u = V + Us, (D9) 

so do the density fluctuation and bond angle variables: 

9 = 60 + 9^; 6p = 6po + 6ps. (DIO) 

The defect source term, S'(r), can be expressed in terms of ps by noting first that, for 
a distribution of dislocations with Burger vector density b{f) = J2ab°'^{^ ~ ^a), S{r) = 
—V X b{r) from Eq. (p4]) , and hence, 

V-'Sir) = -i-V ■ Us = (Dll) 
i — a i — a 

where a = A/(2/i + A) is the Poisson ratio. 

The partition function of annealed dislocations |Q involves averaging over thermally 
excited Burger's vectors with weight proportional to 
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exp(^-E,J d\b{r) ■ 6(f)) = exp(^-E, J rfV(4(V^,)2 + (^^^(V5p.)')) • (D12) 
Thus, from Eqs. ( |D8| - |I7T2D , it follows that 

^annealed = jv5 PqV5 psVOoVe s exp ^ _ "^""y "^^'^ ^ ^13) 



with 

J^annealed= ^ j (fri^AjlOl + (2/i + A)(5po)^) + / d^r{6 po)i6 Ps) 

+ 2(1^ / dMSpsf + E^I dh-(Amr + j^l-y^{V5p.yy (D14) 

Consider now the linear response of the system to a probe coupled to the bond angle. 
One finds easily that 



p{q) p E^q^ 

showing that the shear modulus p vanishes at long wave lengths for finite E^. (i.e., above 
the melting temperature). Similar behavior for p{q) was found by Marchetti and Nelson in 
a dislocation loop model of the melted Abrikosov flux lattice [^. In contrast, the coupling 



between 5po and 5ps prevents the bulk modulus from vanishing at long wave lengths. 
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FIGURES 

FIG. 1. (a): Directed lines in 1+1 dimensions subjected to random point pinning due to a 

disordered substrate, (b): Tethered network of particles in 2+0 dimensions subjected to random 
point pinning due to a disordered substrate. 

FIG. 2. (a): Annealed dislocation disorder embedded in an otherwise six-fold coordinated 
membrane. Heavy lines join the 5- and 7-coordinated sites at the cores of the dislocations. The 
disordered substrate potential is not shown. 

(b): Random substitutional disorder in a polymerized membrane which preserves the 
sixfold coordination of a perfect lattice. The disordered substrate potential which acts on 
this lattice is not shown. 
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